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Abstract 

We determine the qq condensate for quark masses from zero up to that of the strange 
quark within a phenomenologically successful modelling of continuum QCD by solving the quark 
Schwinger-Dyson equation. The existence of multiple solutions to this equation is the key to an 
accurate and reliable extraction of this condensate using the operator product expansion. We 
explain why alternative definitions fail to give the physical condensate. 
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1 Introduction 



Low energy hadron dynamics is controlled by the breaking of chiral symmetry, directly 
reflecting the non-trivial structure of the QCD vacuum. This vacuum is dominated by 
long range correlations between quarks and antiquarks. It is the scale of this qq condensate 
that determines the mass of constituent quarks and hence the masses of all light hadrons. 
Recently experiments involving tttt interactions at and close to threshold have allowed the 
uu, dd condensates to be extracted and these confirm a size of ~ —(235 MeV) 3 anticipated 
from phenomenology pQ. How does this value depend on the current mass of the quarks 
within QCD? 

The Schwinger-Dyson equations are the natural tool for investigating strong coupling 
QCD in the continuum. To study the momentum dependence of the quark propagator one 
needs to know the product of the quark-gluon vertex and the gluon propagator. Models 
for these have been developed over the last decade with remarkable phenomenological 
success. The models have sufficient interaction strength below 1 GeV that the vacuum 
becomes non-trivial and a mass gap is created [21 [3]. Consequently, quarks with zero 
current mass have a dynamically generated mass. Indeed, the momentum dependence of 
the quark propagator for any current mass is determined. The aim of this paper is to 
extract from this dependence the behaviour of the value of the qq condensate beyond the 
limit of vanishing quark mass. 

The interest in the value of such a condensate arises in the context of QCD sum-rules. 
There the Operator Product Expansion (OPE) is used to approximate the short distance 
behaviour of QCD. In studying currents like that of ^7^(75)^, with qi = s and qj = u, d, 
the vacuum expectation values of uu, dd and ss operators naturally arise [HE], [6]. The value 
of the qq condensate for the u and d quarks now well determined to be —(235 ± 15 MeV) 3 
by experiment [7j is the result in the chiral limit. However, in the OPE it is the value of 
the condensates away from this limit that actually enters. Since the current masses of the 
u and d quarks are only a few MeV, the resulting condensate is expected to be close to 
its value in the chiral limit, but how close? For the first 20 years of QCD sum-rules their 
accuracy was never sufficient for it to matter whether this difference was even a 10% effect. 
This equally applies to the estimate by Shifman, Vainshtein and Zakharov [H [9] that the 
ss condensate was (0.8 ± 0.3) of the uu and dd values. It is the greater precision brought 
about by the studies of Refs. [U[l0j[TT], for instance, that motivate the need to learn about 
how the qq condensate depends on the current quark mass. Our aim here is to illustrate a 
method for determining this dependence. 

For light quarks, u, d and s, studying the Schwinger-Dyson equation for the fermion 
propagator in the continuum is essential (see for instance p2]), until computation with 
large lattice volumes with small spacing become competitive. 



2 



2 Schwinger-Dyson Equations 



Our aim is to calculate the mass function of the quark propagator for a range of current 
masses. The starting point is the renormalized Schwinger-Dyson equation for the quark 
propagator as depicted in Fig. [TJ 

Sf\p) = Z 2 

x^S F (k)T u (k,p)D tMU (p-k) . (1) 

In the Landau gauge we note that Z\ — 1. The inverse propagator Sp 1 (p) is specified by 
two scalar functions A and Ai: 

S~\p) = A{p 2 ){p l + M{p 2 )) . (2) 

While A is also a function of the renormalisation point \x and so strictly A(p 2 ,fi 2 ), the 
quark mass function M(p 2 ) is renormalisation group invariant. Of course, Eq. ([1]) involves 
not just the quark propagator, Sp, but the full quark-gluon vertex, T v , and the gluon 
propagator D^ v . Rather than solving for the coupled quark, ghost and gluon system, one 
may employ some suitable ansatz for the coupling and interaction in Eq. ([I]) which has 
sufficient integrated strength in the infrared to achieve dynamical mass generation. There 
have been many suggestions for this in the literature [21 E] which have been extensively 
studied. Following the lead of Maris et al. [21 [5J, we will employ an ansatz for g 2 D fJiU (p — k) 
which has been shown to be consistent with studies of bound state mesons, and consider 
other modellings elsewhere. Since this simple model assumes a rainbow vertex truncation, 
the solutions are not multiplicatively renormalisable and so depend on the chosen renor- 
malisation point. The renormalisation scheme is one of modified momentum subtraction 
at some point fi. This we take to be at 19 GeV to compare with earlier studies [15] . We 
will later evolve \i to the more common 2 GeV scale in the MS'-scheme. We use: 

gJ^(?)-a(g 2 )4°J(g) (3) 

where the coupling is described by [21I5J: 
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Figure 1: Schwinger-Dyson equation for the quark propagator 
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a(q 2 ) = - 6 Dq*eM-q 2 /^) 
+ 



log (r+{l + q 2 /^ Q c D T) 
x [l-exp(-g 2 /[4m t 2 ])] , (4) 



with 



m t = 0.5 GeV , r = e 2 - 1 

7m = 12/(33 -2iV» , A Q cd = 0.234 GeV. 

We work with both the Nf = 0,4 limits. The precise value of Aqcd is irrelevant for our 
current study. The parameters u and D control the infrared behaviour of the coupling 
a(q 2 ) for momenta less than Aqcd, or strictly speaking m t . The pion decay constant 
correlates these parameters. Consistent with this we choose uj = 0.4 GeV, D = 0.933 
GeV 2 as a typical parameter set. Solutions are obtained by solving the coupled system of 
fermion equations for A and M. of Eq. (J2J), which we may write symbolically as: 

A{p\ii) = Z 2 (j/,A)-£jj(p,A) , 

(5) 

M(p 2 )A(p 2 ,n) = Z 2 (fji,A)Z m m R (n) + E s (p,A) . 

The £5 and correspond to the scalar and spinor projections of the integral in Eq. (1T1). 
For massive quarks we obtain the solution M (later called M + ) by eliminating the renor- 
malisation factors Z 2 , Z m via: 

Z 2 (fji,A) = l + S D (/i,A) , 

(6) 

7 (., \\ _ _A (^ A ) 

J Z 2 ( M ,A) Z 2 0i l A)m fl ( M )" 

The momentum dependence for different values of ttir are shown in Fig. [2J Our purpose 
is to define the value of the qq condensate for each of these. 

At very large momenta the tail of the mass function is described by the operator product 
expansion so that 



M(p 2 ) asym = m [log(p 2 /A 2 ]) 
2tt 2 C 



+ 



3 p 



2 



\\og{p 2 /Al) 



(7) 



where the first term corresponds to the explicit mass in the Lagrangian, so m is related 
to the quantity rriR(fi) via some renormalisation factors. The second term gives the lowest 
dimension vacuum condensate, where C = (qq). This provides an excellent representation 
of all our solutions. If we included the expression to all orders then the scales Ai and A2 



4 




p 2 (GeV 2 ) 



Figure 2: Euclidean mass functions for different current masses, specified at fi = 19 GeV 
as labelled. The plot illustrates how on a log-log plot the behaviour dramatically changes 
between a current mass of and 3 MeV. These results are essentially the same as found 
by Maris and Roberts [3]. 



would both be equal to Aq C d- However, the leading order forms in Eq. ((Tj) absorb different 
higher order contributions into the two terms and so A x and A 2 are in practice different, 
as we will discuss below. For large masses the condensate piece, C, is irrelevant and so it 
is the leading term that describes the mass function well. In contrast in the chiral limit, 
when m = 0, the second term of the OPE describes the behaviour of the mass function. 
This then accurately determines the scale A 2 . Indeed, its value is equal to IVqcd- We can 
then in the chiral limit easily extract the renormalisation point independent condensate, 
C = (qq), from the asymptotics. 

For non-zero current masses, one can attempt to fit both terms of the OPE in Eq. (jTJ) 
to the tail of the mass function, M of Fig. [2j Comparing the full mass with m q ^ with 
that in the chiral limit, one sees how very small is the contribution that any condensate 
makes to the tail. So while a value for the condensate can be extracted, this procedure is 
not at all reliable. This is because of the difficulty in resolving the two functions in the 
OPE from one another and in fixing the appropriate scales, Ai and A 2 . 
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Strictly in the chiral limit, we may also extract the condensate using: 

-(qq) u = Z 2 (ji,A)Z m (ji,A)N c tT D [ ^Ls F (k,n), (8) 



(2nY 

with (qq) the renormalisation point dependent quark condensate. At one-loop, this is 
related to the renormalisation point independent quark condensate: 



<99>«= UtogTT— ■ ( 9 ) 




which we compare with the asymptotic extraction to good agreement. 

However, for non-zero quark masses, even for very small masses where the condensate is 
still expected to play a sizeable role, we cannot apply Eq. (jSJ) , since it acquires a quadratic 
divergence, cf. Eq. (|7|). Indeed, it is the elimination of this divergence that inspired the 
definition proposed by Chang et al. [T3], which, as we will mention below and explicitly 
show elsewhere, is unfortunately not equal to the condensate of the physical mass function 
and is therefore ambiguous. Consequently, we need a different definition, one close to the 
OPE, Eq. 0. 

This comes about by recalling that the Schwinger-Dyson equation, Eq. (JTJ, has multiple 
solutions. In the chiral limit, there exist three solutions for Sp(p) and its mass function 
M(p 2 ). These correspond to the Wigner mode (the only solution accessible to perturbation 
theory), and two non-perturbative solutions of equal magnitude generated by the dynamical 
breaking of chiral symmetry. These we denote by Sp' w (p), where: 

M w (p 2 ) = 

M(p 2 ) = { . (10) 

M ± {p 2 ) = ±M°{p 2 ) 

Analogous solutions to these exist as we move away from the chiral limit. The existence 
of these is essential to our ability to determine the qq condensate beyond the chiral limit, 
up to the strange quark mass. While here we will present the results for the Maris- Tandy 
model, elsewhere we will provide the details for other more sophisticated modellings of 
QCD and for the extrapolation of lattice results to the continuum, which within the larger 
uncertainties inherent in these procedures confirm the results presented here. 

Crucially, for a given tur we obtain the M~ and M w solutions by inserting the same 
Z2 and Z m found for the M + solution. This ensures that differences in the dynamics of 
the three systems do not influence the ultraviolet running of the current-quark mass in the 
context of the subtractive renormalisation scheme used by Maris and Tandy. The iteration 
process is performed using Newton's method. For the solution M w this is mandatory, since 
it corresponds to a local maximum of the effective action and is therefore not accessible 
using the conventional fixed point iteration scheme. We renormalise at p, — 19 GeV to 
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Figure 3: Momentum dependence of the three solutions M ± (p 2 ) and M w (p 2 ) for a quark 
mass m(/i)=16 MeV, fi = 19 GeV. 

compare with the work of Ref. [13]. However, in the plots that follow we show the quark 
condensates as functions of the quark mass at \i = 2 GeV. This is obtained by one loop 
running from that at 19 GeV and just provides a convenient mass to plot. Later we will 
transform these momentum subtraction results to the commonly used MS scheme at 2 
GeV. A representative example of the solutions is shown in Fig. [31 

As we illustrate below the solutions M~ and M w only exist below some critical mass 
m cr at which a bifurcation occurs. For zero active flavours and u = 0.4 GeV, we obtain 
m cr = 43.4 MeV at [i = 2 GeV. The exact value of m cr is found to be model-dependent. Of 
course, the physical solution M + exists for all quark masses as shown in Fig. [2j If one takes 
the difference between any pair of solutions of M^, M w , then in each of these the first 
term in the OPE, Eq.([7j), explicitly cancels and each has a behaviour qualitatively like M + 
in the chiral limit, Fig. [2} being dominated by a condensate term. We can then determine 
this condensate, C in Eq. (J7J), by integrating any of these differences, like Sp — Sp, as in 
Eq. (JED- Though such differences give a convergent integral, its value is just the difference of 
the condensates, for instance C + — C~ and this is not a measure of the physical condensate 
as Chang et al. presume. This is why we do not use a convergent integral to determine 
the physical condensate, but rather the OPE as we now describe. 
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Figure 4: Condensate extracted through simultaneous fitting of the three solutions to 
the fermion mass-function in the Maris- Tandy model with Nf = and to = 0.4 GeV as 
functions of the current quark mass evolved to 2 GeV in a momentum subtraction scheme. 

3 Extracting the Condensate 



Instead of one single solution, we now have three solutions to the same model, each with 
identical running of the current-quark mass (the first term in Eq. ([7])) in the ultraviolet 
region and differing only by their values of the condensate. Thus, for m^(/i) < m cr , it is 
possible to fit Eq. ([7]) simultaneously to the three mass functions M , M w . The scales A x 
and A2 are found to be the same for any current mass within the given model. A2 is equal 
to Aqcd, while A 1 is roughly twice as big. Its exact value depends on the parameters of 
the model, like uj and D of Eq. (j3J). The condensates and C w are then determined in an 
accurate and stable way. This fitting is performed using a modified Levenberg-Marquardt 
algorithm. The results are shown in Figs. HI [5] for Nf = 0,4. The error bars reflect the 
accuracy with which the mass functions represented by two terms in the OPE expression, 
Eq. (J7J), are separable with the anomalous dimensions specified. 

We see that within errors the condensate is found to increase with quark mass. This 
rise at small masses was anticipated by Novikov et al. [14] combining a perturbative chiral 
expansion with QCD sum-rule arguments. That the chiral logs relevant at very small m q 
are barely seen is due to the quenching of the gluon and the rainbow approximation of 
Eq. ([T]). As we will show elsewhere this has little effect when we model more complex 
interactions including matching with the lattice. 
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Figure 5: Condensate for Maris- Tandy Model with Nf = 4, uj = 0.4 GeV as a function of 
current quark mass evolved to 2 GeV, to be compared with Fig. HI 



We see in Figs. HJ El that the M~ and solutions bifurcate below m cr ~ 43.4(44.0) 
MeV with uo = 0.4 GeV for Nf = 0(4) respectively. But what about the value of the 
condensate for the physical solution M + beyond the region where M~ and M w exist, i.e. 
m#(/i) > m cr l Having accurately determined the scales Ai and A 2 in the OPE of Eq. ([7]) 
in the region where all 3 solutions exist, we could just continue to use the same values in 
fitting the physical M + solution alone and find its condensate. However, this would make 
it difficult to produce realistic errors as the quark mass increases. 

As soon as one allows for solutions for the fermion mass-function that are not positive 
definite, one exposes a whole series of variants on the solutions M~, M w we have already 
considered. Thus there are noded solutions, which have also been discovered recently in 
the context of a simple Yukawa theory by Martin and Llanes-Estrada [T5]. We illustrate 
this within the Maris- Tandy model, for instance with Nf = 4 and uj = 0.4, in Fig. [61 
There the four solutions we have found are displayed. It is interesting to note that this 
noded solution is not limited to the same domain that restricts M~ and M w . These 
noded solutions do develop a singularity in M{p 2 ) beyond m = 51.4 MeV at \i = 2 GeV. 
However, this is compensated for by a zero in A(p 2 ), Eq. (El), until m = 66.3 MeV. Thus 
there exists a solution with a well-defined ultraviolet running of the quark mass exactly 
as the M + solution, as far as m = 66.3 MeV. While at small quark masses we have all 
four solutions, at larger masses there are still two. Consequently, we can accurately fix the 
scales Ai and A 2 of Eq. at each m cr , and so determine the condensates as shown in 
Fig. [71 Indeed, fitting the M + and M^ oded at each value of m#(/i) with common scales in 
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Figure 6: Momentum dependence of the 4 solutions for the fermion mass-function in the 
Maris-Tandy model with m = 20 MeV at ^ = 19 GeV, n F =4, u = 0.4 GeV. 
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Figure 7: Current quark mass dependence of the condensates for Maris-Tandy model with 
tif = 4, uj = 0.4 GeV, including the noded solution of Fig. [6j 
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Figure 8: Condensate for Maris- Tandy Model with Nf = 4, u = 0.4, 0.45, 0.5 GeV as a 
function of current quark mass defined at 2 GeV in MS scheme. 



the OPE equation, Eq. ([7j) allows the condensate for the physical solution to be found for 
much larger quark masses, as shown in Fig. [7J Indeed, these fits confirm that A x and A 2 
are independent of m^(/i). We can then fit the remaining M + solutions shown in Fig. [2] to 
give the physical condensate shown in Fig. [7J for acceptable values of u as determined by 
[T6] . In Fig. [8] we have in fact scaled the quark mass from /i = 2 GeV in the (quark-gluon) 
MOM scheme by one loop running to the MS scheme at 2 GeV using the relationship 
between Amom and A-^g for 4 flavours deduced by Celmaster and Gonsalves [17]. In this 
latter scheme the strange quark mass is ~ 95 MeV |18j . 

Within the range of the Maris- Tandy modelling of strong coupling QCD we find that 
the ratio of the condensates at the strange quark mass to the chiral limit is 

(^) m( M5)=95Mev/(w) m= = ( ^ ± 0-2 ) 3 . (11) 

in a world with 4 independent flavours. In a longer paper we will consider the Schwinger- 
Dyson equations with more sophisticated ansatze incorporating the quark, gluon and ghost 
interactions both in the continuum and on the lattice. Moreover, here all the quarks have 
the same mass and there is no mixing between different hidden flavour pairs. Elsewhere 
we will illustrate the change that occurs in solving the quark Schwinger-Dyson equations 
with 2 flavours of very small mass m Uj d and 1 flavour with variable mass. Of course, in the 
quenched case quark loops decouple and exactly replicate the results given here. Chang et 
al. have proposed that the fact that the M~ and M w solutions only exist for some domain 
of quark masses, m q G [0, m cr ], is directly linked to the domain of convergence of the chiral 
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expansion [191 [20]. While the existence of multiple solutions for the fermion mass function 
is essential for the extraction of the condensate, we note that distinct domains exist for the 
different solutions and that the simplest noded solution exists in a larger domain. Though 
the value of m cr is indeed some measure of the range of validity of the chiral expansion, 
the fact that m cr is both strongly model and solution dependent indicates its value is no 
more than a guide and not likely to be the exact bound claimed by Chang et al. 

What we have shown here is that there is a robust method of determining the value of 
the qq condensate beyond the chiral limit based on the Operator Product Expansion. Of 
course, as the quark mass increases the contribution of the condensate to the behaviour 
of the mass function, Fig. [2J becomes relatively lass important and so the errors on the 
extraction of the physical condensate increases considerably. Nevertheless, the method is 
reliable up to and beyond the strange quark mass. Alternative definitions are not. 
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